
Parameterized Spectral Distributions for Meson 
Production in Proton-Proton Collisions 


John P. Schneider, John W. Norbury, and Francis A. Cucinotta 


(NASA-TM-4t675) PARAMETER I ZEC N95-30337 

SPECTRAL DISTRIBUTIONS FOR MESON 
PRODUCTION IN PROTON-PROTON 

COLLISIONS (NASA. Langley Research Unclas 

Center) 20 p 

hi/73 0057643 


June 1995 




NASA Technical Memorandum 4675 


Parameterized Spectral Distributions for Meson 
Production in Proton-Proton Collisions 


John R Schneider and John W. Norbury 
University of Wisconsin • La Crosse, Wisconsin 

Francis A. Cucinotta 

Langley Research Center • Hampton, Virginia 


June 1995 


Acknowledgment 


John P. Schneider would like to thank R. T. Wheeler of the University 
of Wisconsin for his assistance in implementing Gaussian quadrature. 
This work was supported by NASA grant NAGl-1 154, the University 
of Wisconsin, La Crosse, Wisconsin, and the Wisconsin Space Grant 
Consortium. 


Available electronically at the following URL address: http://techreports.larc.nasa.gov/ltrs/ltrs.htmi 


Printed copies available from the following: 

NASA Center for AeroSpace Information 
800 Elkridge Landing Road 
Linthicum Heights, MD 21090-2934 
(301)621-0390 


National Technical Information Service (NTIS) 
5285 Port Royal Road 
Springfield, VA 22161-2171 
(703) 487-4650 



Abstract 


Accurate semiempirical parameterizations of the energy^differential cross sec- 
tions for charged pion and kaon production from proton-proton collisions are 
presented at energies relevant to cosmic rays. The parameterizations, which depend 
on both the outgoing meson parallel momentum and the incident proton kinetic 
energy, are able to be reduced to very simple analytical formulas suitable for cosmic- 
ray transport through spacecraft walls, interstellar space, the atmosphere, and 
meteorites. 


Introduction 

In order to solve cosmic-ray transport problems, a 
detailed understanding of the proton-proton interaction is 
needed. Galactic cosmic rays consist of approximately 
85-percent protons (ref. 1) and solar cosmic rays consist 
almost entirely of protons and a-particles, whereas geo- 
magnetically trapped particles in the Van Allen belts 
consist of protons and electrons (ref. 2). The interstellar 
medium itself consists primarily of hydrogen (protons) 
and helium (ref. 3). Hence, the proton-proton collision is 
the most common and most numerous interaction 
between cosmic rays and the interstellar medium, as well 
as within the solar system and the Earth’s atmosphere. 

Proton transport can be adequately modeled by using 
the Boltzmann equation in the straight-ahead approxima- 
tion (ref. 2), which is given as 

+o]<t>U.£) = f(E,E')<^{x,E')dE' (1) 

where S(E) is the proton stopping power, a is the media 
macroscopic cross section, (|)(;c, E) is the particle flux 
density, J{E, E') is the secondary-particle differential 
cross section in terms of secondary-particle energy (£) 
and proton energy {E), and jc is a scaling variable. This 
straight-ahead approach is used in current transport 
codes for particle propagation through various media 
(refs. 4-7). 

One major shortcoming of the work done to date is 
that although the secondary-particle-production cross 
section E) in equation (1) is taken to include all pos- 
sible secondary particles, codes have not yet been run 
using meson cross sections (ref. 2). In this paper, param- 
eterizations of cross sections for charged pions and kaons 
suitable for use in equation (1) will be provided. 

The numerical approach currently used in solving 
equation (1) implements the cumulative energy spec- 
trum, which is related to f^E, E) above as follows 
(ref. 2): 

F(r,r') = j^'f(E,E')dE (2) 


where the residual range (r) is given by r= J dElS(E). 
For a very thorough analysis of proton transport and a 
detailed explanation of the numerical procedure for solv- 
ing equation (1), the reader is referred to Wilson etal. 
(ref. 2). Such an approach requires that the spectral distri- 
bution fiE, E) be calculated a large number of times, 
thus making it impractical to utilize complicated analyti- 
cal formulas for f{E, E). Hence, the parameterizations 
derived herein will be presented in the simplest possible 
terms (without losing accuracy, of course). 

An interaction cross section can be written in various 
forms, and high-energy physicists are interested in the 
Lorentz-invariant form E d^ofd^p (where a is the total 
cross section and p is an energy momentum), which can 
be calculated from first principles using quantum field 
theory and is thus a convenient way to check theory with 
experiment. In this light, all early theoretical and experi- 
mental work done with mesons is presented in Lorentz- 
invariant form (refs. 8-11). 

However, f{E, E) from equation (1) is the energy- 
differential cross section do/dE (or spectral distribution). 
The Lorentz-invariant cross section can still be used as it 
contains all the information one might need (including 
angular and spectral distributions as well as the total 
cross section). It now becomes a matter of retrieving this 
information. 

In one of the early papers (ref. 1 1 ), all then-available 
accelerator data were gathered, and a very accurate 
parameterization of the Lorentz-invariant cross section 
was fit to the data. The aim of the current work is to uti- 
lize this representation of the Lorentz-invariant cross sec- 
tion in order to generate a representation of the spectral 
distribution that will then be parameterized as a simple 
function of incident proton energy and secondary meson 
momentum. The final result will provide the proper dif- 
ferential cross section for equation (1). 

Symbols 

A I parameter from table 1, mb/(GeV /c ) 

^ 2 , C, C 2 parameters from table 1, (GeV/c)“^ 

C[, y parameters from table 1, dimensionless 



C3 parameters from table 1 , (GeV/c)“^ 

c speed of light 

c.m. center-of-mass frame 

d^O differential cross-section element, mb 

d p differential momentum element, GeV/c 

E energy, GeV 

/ continuous spectrum /(£, E) from equation (2) 

K kaon 

M = reaction, A ^ B ^ C + D), 

GeV/c^ 

m mass, GeV/c 

Pj^ energy-momentum four vector 

p momentum vector, GeV/c 

q parameter from equation (5) 

r residual range from equation (2) 

S{E) proton stopping power from equation (1) 

s Mandelstam energy variable 

X scaling variable used in representation of 

Lorentz-invariant cross section 

X{* ratio of c.m. momentum to maximum 

transferable momentum 

a, P, X parameters in equation (14) 

0 scattering angle 

a cross section, mb 

(|) particle flux density from equation (1) 


s is the Mandelstam energy of the system, and is the 
secondary-particle rest mass. 

In terms of this variable, the following parameteriza- 
tions are then given (ref. 1 1) for pions: 


E 

d^p 


— - (1 -x)’exp 

1 + (4mp/i) 


-^2Pl 
1 + (4mp/s) 


where 


(4) 


C] + C2P_i + C^p\ 

1 + i4mps) 


and for kaons: 


(5) 


3 

( 6 ) 

d^p 

where the parameters A|, A2, C, Cj, C2, C3, and y are 
given in table 1 , and nip is the proton rest mass. 

Now, if the spectral distribution is to be written in 
simple terms (i.e., constants, mass, momenta, and ener- 
gies of the interacting particles), then it will be desirable 
to first express the Lorentz-invariant cross section in 
such terms. Specifically, the Mandelstam energy variable 
(s) and the fractional momentum (xp ) must be simply 
represented. 

The Mandelstam representation of the energy for any 
system A+B^C-^D is defined as 


Subscripts: 
max maximum 

p proton 

s secondary particle 

7C pion 

± perpendicular 

// parallel 


(7) 

where and are energy-momentum four vectors, 
given by Pf^ = Ef^lcpp^, where the bold-faced p^ repre- 
sents the usual three-dimensional momentum vector. 
This can then be rewritten as (see appendix A) 

s - {lE^ntglc^) (8) 


Invariant Cross Section 

The Lorentz-invariant cross section is given in terms 
of the following scaling variable (refs. 1 1 and 12): 


For the current problem, A and B are protons and is 
the total incident proton energy (given by the sum of the 
proton kinetic energy (T|ab) mass (nip). Thus, 

the following representation is given: 


X 



(pi+"«5) 


1/2 


(3) 


= 4/n^ + 2m^r,^b 


(9) 


where Xp* is the ratio of the parallel component of the 
center-of-mass (c.m.) momentum to the maximum trans- 
ferable momentum, p^_ and pn are the perpendicular and 
parallel components of the c.m. momentum, respectively, 


The fractional momentum x// is given by p// 
Pc.m.,max ^ where Pc.m.,max the maximum trans- 
ferable momentum allowed by the kinematics in the 
center-of-mass frame. For the arbitrary system 
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Table 1. Parameters From Equations (4), (5), and ( 6 ) Used in Representation of Invariant Cross Section 


[Data based on ref. 1 1] 


Particle 

A], mb/(GeV 2 /c 3 ) 

Aj, (GeV/c)’' 

Y 

C, (GeV/cr' 

C, 

Cj, (GeV/c)“‘ 

C 3 , (GeV/c)“^ 

71+ 

153 

5.55 

1 


5.3667 

-3.5 

0.8334 


127 

5.3 

3 


7.0334 

^.5 

1.667 


8.85 

4.05 


2.5 




K" 

9.3 

3.8 


8.3 





A+5->C + D, Pc.m.,max given by (ref. 13 and 
appendix A) 


^c.m.,max 


, 2 4 ..2 4.2 2,. 2 8 ^ 1/2 

(s - nif^c -Me) -Am^M c 

As 


( 10 ) 


where 




( 11 ) 


Once again, for the current problem, this simplifies to 


c.m.,max 


, 2 4 , 2 4.2 ,, 2 2 8 -, 1/2 

{s-m^c -4ntpC ) - Ibrn^nipC 


As 


( 12 ) 


Hence, this allows the Lorentz-invariant cross sec- 
tions from equations (4) and (6) to be expressed solely in 
terms of the masses, momenta, and energies of the parti- 
cles (and the parameters given in table 1). 


Spectral Distribution 

The spectral distribution doldE is expressed in terms 
of the invariant cross section via (ref. 1 1 and appendix A) 

^ = - f ^4^ dipl) (13) 

dE Pn j ^3^ 

By inserting the invariant cross sections from 
equations (4) and (6) into equation (13) and integrating, a 
representation of dddE will be produced in terms of pro- 
ton kinetic energy (Ti^^) and meson parallel momentum 
(p//). Because of the nonlinear exponential form of 
equations (4) and (6) and the nonlinear integration in 
equation (13), the integrals cannot be performed analyti- 
cally. A numerical procedure utilizing adaptive Gaussian 
32-point quadrature was implemented (see the computer 
program in appendix B); this procedure produced the 


data points seen in figure 1, which shows the spectral 
distribution for both positively and negatively charged 
pions and kaons at a variety of proton energies. 

By fitting a function of pn and to the data in 
figure 1, the parameterization will be solved. Clearly, the 
data are well represented by decreasing the exponentials, 
and the solid-line curves in the figures represent the best 
fits obtained by utilizing the following exponential func- 
tion to represent doldE: 

^ = —ac\p{-^pff) (14) 

dE Pn 

where a, P, and % vary with respect to the incident proton 
kinetic energy (T|^i,) and must be parameterized as such. 
This procedure is accomplished by determining which 
values of a, p, and x best match the various proton ener- 
gies from equation (14) with the output from 
equation (13), and then by finding functions that accu- 
rately provide the same values if given simply the proton 
kinetic energy (Tia^). The following functions were 
found to be extremely accurate: 

For a. 


a = a, In(a2 + 7']3b) +“3 1" (“4 + ^lab) ^ + “5 (1^) 


forp. 


p = p,r^ ln(p3 + 7',,b)“^'-P57'i (16) 

and for X, 


7fab + Z3 
X = X\ 

T^b + Xs 


(17) 


where the subscripted parameters of equations (15), (16), 
and (17) are given in table 2. The accuracy of 
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Table 2. Parameters for Representations of a, P, and X Used in 
Equations (15), (16), and (17) 


Parameter 

7t^ 

n 


K" 



1.6209 

2.0555 

0.513 

1.0023 


«2 

-.0992 

2.5949 

3.1993 

15.612 

a 

«3 

-.0958 

-.546 

-.0421 

-.0711 


«4 

-.0217 

689.37 

-.8582 

2.1746 


^5 

2.5884 

21.24 

-.6736 

-2.7214 


Pi 

57.4728 

72.1723 

25.4335 

66.5962 


P2 

-.6054 

-.6595 

-.5313 

-.5619 

p 

P3 

.7049 

.7049 

-.1231 

-.122 


P4 

-.534 

-.5335 

-.9639 

-.9198 


Ps 

42.169 

51.658 

12.887 

28.343 


Pa 

-.7979 

-.8771 

-.8233 

-.8742 


Xi 

1.9096 

1.8773 

2.3452 

33.621 


X2 

1.6939 

1.4223 

1.1829 

-.0301 

X 

1 

.5613 

.4237 

.5525 

-.5638 


X4 

1.7668 

1.4898 

1.2645 

-.5344 


X5 

.6666 

.5082 

.7275 

6.361 


equation (14) is evidenced by the very fine agreement 
between the data points and the solid-line curves seen in 
figure 1. 

Results 

The agreement between the fitted approximation 
using equation (14) and the output from the numerical 
quadrature used to evaluate the integral of equation (13) 
is found to be extremely accurate for a wide range of pro- 
ton energies relevant to cosmic-ray interactions. 

A final interesting point in regard to Nagamiya and 
Gyulassy’s representation (ref. 13) of the maximum 
transferable momentum (eqs. (10) and (12)) is that the 
representation involves a square root that allows a lower 
bound to be set on ^lab- below which the root becomes 
imaginary. Because the only variables besides Tiat, that 
are involved in the square root are the masses of the par- 
ticles (ref 14), the following lower limits on are 
obtained for meson production in proton-proton colli- 
sions. Thus, for pions, 

Tiab > 0.29 GeV 


and for kaons, 

7'i,b>1.12GeV 

which are just the threshold energies. 

Concluding Remarks 

A parameterization of the energy-differential cross 
sections for charged pion and kaon production from 
proton -proton collisions has been derived that accurately 
reproduces results obtained by numerically integrating 
the Lorentz-invariant cross section. Because the parame- 
terization is simply a function of incident proton kinetic 
energy and outgoing meson parallel momentum, the 
results can be easily reproduced and should be very use- 
ful in computer transport codes that have not yet been 
tested with meson cross sections. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
March 30, 1995 
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Appendix A 

which yields 


Proofs of Formulas Used in Text 

Ec^Ejy = 

When relating enerev and momentum via 

(Al) 

Mandelstam Energy Vanable 

The Mandelstam energy variable {s) for the reaction 



A -h B — ^ a ^ D 

eI = i-q)^c^ + mlc^ 


by definition, is given as ^ 


the difference between E^ and E^ is given as 


(Pa 

eI-eI = 


In lab frame, Pr-0 and En = rrioc^. Thus, s is 

expressed as 

Now, since 


{Pa +Pb)^^^^ 



= (Pa Pb'^ ^PaPb) 

we have 


= + 2p^Pg) /c^ 

^C~^D ~ ('^c ” ^d) 


~ ^A ^B ^ (2/c ) P aP B 

where + Ep) is given by equation (A1 ). Hence, 

= m^+m^+ (2/c^) {E^/c,Pj^) ■ (Eg/c,pg) 

Ec-Ed = irn}~-m]))c‘^ls^''^ 

(A2) 

= ml+ml+ (2/c2) {E^/c, p^) • (nig/c, 0) 

Spectral Distribution 

The spectral distribution is given as 


= +m^+ {lE^mgfc ) 


d<5 71 f £ d^o , . 2n 


Maximum Transferable Momentum 

' .v» ^ .V 


The maximum transferable momentum for the 

For proof, we let 


reaction 

A. B — ^ C D 

Ed^a _ E d^a 

(A3) 

in the center-of-mass frame where p =Pa~ ~Pb 

d^p <^Pn Pi dp_i dQ 

q = Pc = ~Pd^ is given by 

Integrating with respect to 0 gives 


w = [ (» - '«y - 1 

^ E d^a E d^a 

2K = ; 

(A4) 

where 

d^p ^Pn Pi ^Pi 

M = 

We can simplify further via 


Here, s is defined as 

= Pi + pI 


(PA+PB)^f^^ = 

E=[(pl^pl)c^^my-]^^ 


= [{Ectc,q) + {Ejylc,-q)\^lc^ 

dE _ Plf^ 
dpn E 


= [ + (ffl/c, 0) ] 2/c2 

, E dE 



^Pn = , 

(A5) 

= (£c + ^o)^ 

P//^ 
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and 


d{p\) = '2-P^dp^ 

P±dP± = d {p\) n (A6) 

Substituting equations (A5) and (A6) into the denomina- 
tor on the left-hand side of equation (A4) gives 

(P’a _ E (f-g 

i3 E dE r , . 2 N T 
dp -[dip\)H] 

Ptr 

and 


+ (^C ” ^£)) 

~ ^ ('”c “"*£>) 

= I [5 + /4jc^ 

= [i ^ ^ me - ^ me - m^ - 4meC^] /4jc^ 

= [5^ - 2s{^m^ + + ^me-»i^j^c*]/4^c^ 

Note that q is the transferred momentum and that for 
given particles A, B, and C and given incident energy 
^beam’ variable left in q is m/), which can be no 

larger than + nig. Hence, the maximum transferable 
momentum is given as 


71 E d^a 
c^Pn d^p 


d^g 
dE d 


(pi) 


(A7) 


Finally, integrating with respect to yields the desired 
result 



(A8) 


which allows us to isolate E(^, Ed, and q in terms of 
niD, and s as 


«nu« = {s^-2sc^[ml+ + 

+ [mi - (m^+TOg)^] c* } I2cs^^^ 

and if we let M = /n^ -f- mg, the result is 

^max = 2sc^ {ml + M^) + (mc-Af^)^c^] /2cs^^^ 

= [[ j ■- mlc"^ - - 4mcM^c^] 

which is the desired result. 
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Appendix B 


invcs 


Computer Program of Spectral 
Parameterizations 

The following computer program utilizes the param- 
eterization given by Badhwar, Golden, and Stephens, 
(ref. 1 1 ) to plot a graph of the spectral distribution. The 
integration is performed by using adaptive Gaussian 
32-point quadrature to find the integral. 


mkaon 

mp 

mpion 

particle 

pcmmax 


The program is written in FORTRAN-77, which 
does not support the many fonts found in this paper; 
hence, naming the variables can become quite tedious. A 
complete list of the meaning of all the variables found in 
the main procedure is given as follows: 

A, B, C, Cl, parameterization constants (Aj, A2, C, Cj, 
C2, C3 C2, and C3, respectively) from table 1 

dppar dummy variable used to vary pparcm 

inside loop 

dsdE = doldE, energy -differential cross section 

e total energy of incoming proton 

fac represents one part of equation ( 1 ) that 

does not vary with respect to p/f or 

intgrl value of actual definite integral 


ppercm 


pparcm 


q 

R 


s 

Tlab 

X 


xparsq 

xtilde 


value of invariant cross section, returned 
by function of same name 

mass of kaon 

mass of proton 

mass of pion 

determines which particle the program 
will investigate 

maximum transferable momentum 

= pj^, perpendicular component of c.m. 
momentum 

= p//, parallel component of c.m. 
momentum 

= q, from equation (5) 

= y, parameterization constant from 
table 1 

Mandelstam energy 

kinetic energy of incident proton 

= 1 - xtilde, as seen in equations (3) 
and (4) 

= from equation (3) 

scaling variable found in equation (3) 
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IMPLICIT DOUBLE PRECISION (a-z) 

INTEGER particle 

COMMON mass, A«B,C,Cl,C2«C3«R,q« s,parcicle« £ac,xparsq 

mass ot proton in GeV/c**2 

mp - .93827231 

mass of pion in Gev/c**2 

o«>ion - .1395679 

naaa of kaon in GeV/c**2 

mkaon ■ .493646 

pi * 3.141592653569793238 

do 2000 Tlab « 10. 
do 1000 particle • 1,4«1 

IF (par tide. EQ.l) THEN 
for positive pion. . . 

K - 153. 

B - 5.55 
R -1. 

Cl * 5.3667 
C2 - -3.5 
C3 » 0.8334 
mass » nf>ion 

open ( 50 « f ile«' intppOOlO .dat ' ) 

ELSEIF (particle .eq. 2) then 
for neqative pion. . . 

A - 127. 

B « 5.3 
R « 3 . 

Cl » 7.0334 
C2 - -4.5 
C3 > 1.667 
Biass « npion 

open (50, file*' intnp0010.dat' ) 

ELSEIF (particle . eq . 3 ) then 
for positive kaon. . . 

A • 8.85 
B * 4.05 
C - 2.5 
nass ■ nkaon 

open (50, file* ' intpkOOlO .dat ' ) 

ELSEIF (par tide. eq. 4) THEN 
for negative kaon. . . 

A * 9.3 
B - 3.8 
C * 8.3 
mass * mkaon 

open (50, file*' intnkOOlO .dat' ) 
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ENDIP 


* Mandelstam energy 
e a Tlab mp 

s * 2.*(ni)**2) 2.*e*nV 

• nytiriimim transferable momentum 

pcanax * DSQRT( i (s - mass**2 - 4 . • (mp* *2) ) • *2 - 16 . * (mass**2 ) • (b«)** 2) ) / (4 . *s) ) 
write(6, *)particle,Tlab,e,s« ' pcmmax « \pcasnax 

♦ (neither perpendicular nor parallel conponent can be greater than pcxnmax.) 

• (lower limit on integral is 0. Upper limit is pcmniax**2, since we're 

• integrating with respect to the square of the perpendicular momentum, 

* and changed variables) 

* 'fac' is a COMMON variable, and is used in invcs: 

fac - 1. ♦ 4.*(irp**2)/s 

dppar ■ pcmnax/1000 . 
pparcm « 0. 

do 10 pparcm«0 . ,pcmmax, dppar 

xparsq ■ (pparcm/pcnoax) **2 

intgrl - AdaptiveGauss (OFLOAT(O) ,pcovnax**2) 

IF (intgrl.GE.O.) THEN 

dsd£ « (pi/pparcm) * intgrl 
write (50, *) pparcm, dsd£ 

ELSE 

wite(6, *) 'did not work: ', pparcm, intgrl 
ENDIF 

10 continue 

close (50) 

1000 continue 
2000 continue 

stop 

end 


Adaptive Gaussian Quadrature... 

sent in: lower and upper integration limits (ppercm limits) 

returns the integral approximation 'Approx' of integral 
of Lorentz invariant cross section (invcs) . 

DOUBLE PRECISION FUNCTION Adapt iveGauss (Qllim, Qulim) 

IMPLICIT DOUBLE PRECISION (a-z) 
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HTTEGER z, particle 

COMMON mass« A«B« C« Cl« C2« C3 , R. q, s,particle« fac, xparsq 

DIMENSION low(lOOO) 

DIMENSION mid(lOOO) 
dimension uppUOOO) 

DIMENSION sum (1000) 

DIMENSION tol(lOOO) 

DIMENSION sav(S) 

Approx ■ 0. 

z ■ 1 

tol(z) ■ l./(10.**10) 

low(z) • Him 

mid(z) ■ (ilia uiia)/2. 

upp(z) ■ ulim 

sum(z) « Gauss32 (low(z) , upp(z) ) 

100 IF (z.GT.O) THEN 

si ■ Gauss32 (low(z) «mid(z) ) 
s2 « Gauss32 (Aid(z) , upp (z) ) 

sav(l) « low(z) 
sav(2) « niid(z) 
sav(3) a> upp(z) 
sav(4) « tol(z) 
sav(5) » suffi(z) 

z ■ z - 1 

IF (DABS(sl S2 - sav(5) ) .LT.sav(4) ) THEN 
Approx ■ Approx ^ (si '•> s2) 

ELSE 

IP (z.GE.999) THEN 

AdaptiveCauss « -1. 
goto 9999 

ELSE 

calculate for right-hand side subinterval 
z » z ^ X 
low(z) « sav(2) 
upp ( z ) « sav ( 3 ) 
mid(z) > (low(z) upp(z))/2. 
tol(z) * sav(4)/2. 
sum(z) * Gauss32 ( low( z) « upp< z) ) 
calculate for left-hand side subinterval: 
z » z ♦ 1 
low ( z ) > sav ( 1 ) 
upp(z) » sav(2) 
mid(z) « (low(z) > upp(z))/2. 
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tol(z) * tol(z-l) 

sum(z) Gauss32 ( low ( z ) , upp ( z ) } 

Etn^IF 
ENDIF 
GOTO 100 
ENDIF 

AxiaptlveGauss ■ Approx 

9999 return 
end 


* 32-point Gaueeian Quadrature, used in conjunction w/ AdaptiveQuad above 

DOUBU PRECISION FUNCTION Gauss32 ( Ilia, ulim) 

Is^licit X>oxjble Precision (a-z) 

Integer i, degree 

Dimension WgC(32), Zero(32) 

degree * 32 

CAXX LEGEND (Wgt, Zero) 

sum * 0.0 

DO 300 i«I, degree 

X » (Zero(i) * <ulim - Him) ^ Him > ulim)/2. 

NewWgt * Wgt(i)*(ulim - Him)/2. 

sua » sua > invcs (x) *NewWgt 
300 CONTINUE 

Gauss32 - sum 

RETURN 

END 


Lorcntz-invariant cross section of DSQRT( temp) • • 


DOUBLE PRECISION FUNCTION invcs(teap) 

IMPLICIT DOUBLE PRECISION (a-z) 

INTEGER particle 

COMMON mass , A, B , C , Cl , C2 , C3 , R , q, s , particle , f ac , xparsq 
pperca • DSQRT(temp) 

xtilde - DSQRT( xparsq + 4 (pperca* *2 + mass**2)/s) 

X ■ 1. - xtilde 
if (x.GE.O.) then 

if (particle .LT. 3 ) then 

q a (Cl + C2*ppercm C 3 * (ppercm**2) > /DSQRT( fac) 

invcs - (A/(fac**R))*(x**q)*DEXP(-B*ppercm/fac) 

else 

invcs - A*(x**C) *DEXP(-B*ppercm) 
endif 

else 
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invcs » 0. 
endif 
return 
end 


* Legend: Set initial values for Laguerre Weights and Points.... 

• 

SUBROUTIHS LEGEND(WL,P) 

XiqpliciC Doubl« Pracision (A-H,M,0-Z) 

OiAension WL(96),P(96) 

* L«0u«rre weights are WL, Laguerre points are P 

HL( D- 0.7018610009470012118908133214030-02 
HL( 2)- 0.1627439473090566083440799971S3D-01 

NL( 3)- 0.2S39206530926207136880656811680-01 

HL( 4)- 0. 3427386291302141670817027474530-01 

ViL( 5)« 0. 4283589802222664225800352610920-01 

WL( 6)« 0.5099805926237623281832567201380-01 

WL( 7)- 0. 5868409347853556761592797030860-01 

WL( 8)- 0.6S822222776361790S027S266888SSO-01 

WL( 9>- 0. 7234579410884854795693588158660-01 

WL( 10)- 0.7819389578707024342218989687580-01 

WL( ID- 0.8331192422694676767702670616700-01 
WL< 12)- 0.8765209300440393418163420591330-01 

WL( 13)- 0.9117387869576389775039260143790-01 

WL( 14)- 0.9384439908080452474870014611950-01 

WL( 15)- 0.9563872007927490086209854069920-01 

WL( 16)- 0.9654008851472785021630151902630-01 

WL( 17)- 0.9654008851472781031766157155970-01 

WL( 18)- 0.9563872007927491994405677644410-01 

WL( 19)- 0.9384439900080464270989651254240-01 

HL( 20)- 0.9117387869576393244486212097400-01 

VIL( 21)- 0.8765209300440389948716468637710-01 

WL( 22)- 0.8331192422694674165617456651490-01 

t«L( 23)- 0.7019389570707031801529936387850-01 

WL< 24)- 0.7234579410884842826101603918690-01 

V(L( 25)- 0.6582222277636191713756641519240-01 

NL< 26)- 0.5868409347853561098401486972880-01 

ML( 27)- 0. 5099805926237615735785446702270-01 

HL( 28)- 0.4283589002222673506570949086840-01 

tlL( 29)- 0.3427386291302145313736327025820-01 

WL( 30)- 0.2539206530926200111250579105620-01 

WL( 31)- 0.1627439473090559144546096064300-01 

WL( 32)- 0.7018610009470114250752781348550-02 

P( D- -0.9972638618494816670212443909800^00 

P( 2)- -0.985611511545268423328813867101D-fOO 

P( 3)- -0.964762255587506528664576421761O«00 

P( 4)- -0.934906075937739763914713364557D-.00 

P( 5)- -0.8963211557660521883184046032510-.00 
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P( 

6) 

P( 

7) 

P( 

8) 

P( 

9) 

P( 

10)^ 

P( 


P( 

12)> 

P( 

13) 

P( 

14) 

P( 

15) 

P( 

16) 

P( 

17) 

P( 

18) 

P( 

19) 

P( 

20) 

P( 

21) 

P( 

22) 

P( 

23) 

p( 

24) 

P( 

25) 

p< 

26) 

P( 

27) 

P( 

28) 

P( 

29) 

P( 

30) 

P( 

31) 

p< 

32) 


return 

end 


-0.849367613732570012063938236224D>00 
-0.79448379S9679424410630S8822692D-fOO 
-0.732182118740289711311675091565D+00 
-0 . 6630442669302152433036924839580^00 
-0 . 5877157572407624291921379722230+00 
-0 . 5068999089322294149262404516780+00 
-0 . 42135127613063S3951451S1989027O+00 
-0 . 3318686022821276671024293136720+00 
-0. 2392873622521370403914975355520+00 
-0 . 1444719615827964945276029595790+00 
-0 . 4830766568773830257205181837320-01 
0 . 4830766568773830777622224630360-01 
0.1444719615827964806498151517640+00 
0.2392873622521370438609444875060+00 
0.3318686022821276462857476019510+00 
0.4213512761306353604506824694910+00 
0 . 5068999089322294565596038751210+00 
0 . 5877157572407623459254111253360+00 
0 . 6630442669302152314259046761440+00 
0 . 7321821187402897529450385150080+00 
0.7944837959679424271852710148780+00 
0.8493676137325700120639382362240+00 
0.8963211557660521466850411798080+00 
0 . 9349060759377397639147133645570+00 
0 . 9647622555875064592756373826890+00 
0.9856115115452683678176626358440+00 
0.9972638618494816115100931597230+00 
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Figure 1. 
denote 
curves 



(a) Positive pion production. 



(b) Negative pion production. 

Spectral distribution d^lclE plotted against parallel component of outgoing meson momentum. Data points 
actual values obtained by performing numerical quadrature on integral portion of equation (13); solid-lme 
denote the fit approximation of equation (14). 
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